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. Abstract. The memory function approach to spin dynamics in doped antiferromagnetic 

insulator combined with the assumption of temperature independent static spin correlations 
and constant collective mode damping leads to oj/T scaling in a broad range. The theory 
^ ' involving a nonuniversal scaling parameter is used to analyze recent inelastic neutron scattering 

results for underdoped cuprates. Adopting modified damping function also the emerging central 
peak in low-doped cuprates at low temperatures can be explained within the same framework. 
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Q . 1. Introduction 

It is by now experimentally well established that magnetic static and dynamical properties of 
high-T c cuprates are quite anomalous. Early inelastic neutron scattering (INS) experiments 
on low-doped L^-xSr^CuC^ (LSCO) [U [2] revealed that local, i.e., q-integrated dynamic 
spin response in the normal state (NS) exhibits anomalous oj/T scaling, not reflected in 
instantaneous spin-spin correlation length £y which shows no significant T-dependence below 
room temperature. Subsequently similar behaviour has been found in a number of other 
compounds, i.e., in underdoped YBaCusOg+x (YBCO) and in Zn-doped YBCO [2j. More 
recent INS experiments on heavily underdoped (UD) cuprates, including Li-doped LSCO [3], 
YBCO gl El El E|, and Pr^^LaCe^CuO^^ (PLCCO) [8], confirm the universal features of 
anomalous NS spin dynamics so that oj/T scaling is found in a broad range both in q-integrated 
susceptibility Xl( w ) [3E1IZ] and in Xq( w ) at the commensurate AFM q = Q = (tt, it) [3JE1IH]- 
Typically, Xq( w ) is a Lorentzian with the characteristic relaxation rate scaling as T = aT, 
but with a nonuniversal a [21 El E]- Similarly, Xl( w > T) = x'l^-i ®)f( u /T) has been used [21E1E]j 
with f(x) = 2/-7T arctan[Aix + A2X 3 ] and material dependent A\^- It has been also observed 
that at low T < T g some intensity is gradually transfered into a central peak (CP) [31 0] whereas 
the inelastic response saturates. This freezing mechanism appears to be entirely dynamical in 
origin since £t as well as the integrated intensity are unaffected by the crossover. 

The present authors introduced a theory of spin dynamics in doped AFM |9j which describes 
the scaling behavior as a dynamical phenomenon based on two experimental observations: a) 
£t is (almost) independent of T, and b) the system is metallic with finite spin collective-mode 
damping. Then the system close to AFM naturally exhibits oj/T scaling in a wide energy range, 
with saturation at low-enough T. 

Our starting point in the analysis of recent INS experiments is the dynamical spin 
susceptibilty [9] 



where the spin stiffness r/ q ~ 2 J (~ 240 meV) is only weakly q-dependent (J is the exchange 
coupling), w q = (^q/Xq) 1 ^ 2 is an effective collective mode frequency, x q = Xq( w = 0) is the 
static susceptibility and M q is (the complex) memory function containing information on col- 
lective mode damping 7 q = M''(oS). In the NS of cuprates low-frequency collective modes 



at q ~ Q are generally overdamped so that 7 q > 
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UD cuprates close to the AFM phase 
have low charge-carrier concentration 
but pronounced spin fluctuations whose 
dynamics is quite generally governed by 
the sum rule 
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where C q is strongly peaked at Q with 
a characteristic width = 1/£t- 

Moreover, the total sum rule is for a 
system with local magnetic moments 
(spin 1/2) given by (1/N) £ q C q = (1 - 
c/J/4, where is an effective (hole) 
doping. 

While the formalism so far is very 
general, we now introduce approxi- 
mations specific to UD cuprates 
INS experiments listed above indicate 
that within the NS the effective q 
width of Xq( w )> i- e -, dynamical k{uj), 
is only weakly T- and w-dependent, 
even on entering the regime with the 
CP response [I]. Within further analysis we assume the commensurate AFM response at Q and 
the double-Lorentzian form C q = C/[(q — Q) 2 + Kj,] 2 although qualitative results at low uj do 
not depend on a particular form of C q . 



Figure 1. Relaxation rate Tq vs. T for different 
parameters £ and 7=100 meV. Inset: dependence 
of slope parameter a on £ for two different 7. For 
comparison, an estimate for a is included. 



2. Paramagnetic metal: 

We also assume that the damping 7 q (w) is dominated by particle- hole excitations being only 
weakly q and uj dependent. Hence 7q(w) ~ 7, with 7 a phenomenological parameter. The 
assumption of constant 7 in Eq. (fT]) then leads to 
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Note that recent INS data are fully consistent with this form which has been used to extract 

r Q (T) [5j. 

We next exploit the sum rule, Eq.([2]), to determine T q . As shown elsewhere [10] Tq is 
mainly determined by the parameter 



C = 7T7C Q /(27y), 



(4) 



subject to T <^ 7 which is experimentally relevant. The results are presented in Fig. 1 for a range 
of C = 1 - 8. While r Q (0) = T° Q ~ 7 exp(-2C) [TO] , for T > T° Q the variation is nearly linear 
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Figure 2. Temperature evolution of %q(u;,T) for 
UD YBCO with x = 0.45 [7] (symbols) compared 
with theoretical result, Eq. ([3]), with £ = 1.8 and 
7=60 meV (lines). 



Figure 3. Scaling of theoretical 
normalised T) based on Q and 7 as 

in Fig. ([2]). For comparison, the scaling 
function f(x) as frequently used in fits 
to experiments is also plotted. 
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Tq ~ aT being a manifestation of the 
uj/T scaling. But a is not universal and 
depends on £ (and weakly on 7). Then 
to leading order 

r Q (T) * max{- , 7 exp(-2C)}, (5) 

from which a = 7r/(2£) is identified [10J. 
Note that the above Tq leads to XlO^) 
consistent with the "marginal Fermi 
liquid" model introduced by Varma et 
al. [11 J . However, contrary to the usual 
assumption of proximity to a quantum 
critical point where £t oc 1/T, here £r ~ 
const. 

Recently a number of INS experi- 
ments have been reported where Xo( w ) 
can be well described by Eqj3j includ- 
ing linear-in-T behaviour of Tq where, 
depending on material a ~ 0.18 — 0.75 [3J, El [S] . All these a require rather large £, which implies 
very low saturation Tq, whereas INS data for YBCO and LSCO indicate quite substantial Tq. 
However, saturation of Tq, setting in for T < T g , is accompanied by simultaneous appearance 
of the CP which absorbs the 'missing' sum rule. 

In Fig. 2 INS measurements by Hinkov et al. [6l [7] on UD YBCO with x = 0.45 [7] are 
presented together with theoretical curves were the only relevant parameter is £ = 1.8, while 
in Fig. 3 theoretical scaling function x'q(w, T)/x'q(uj, 0) for the same ( (and 7) is plotted. The 
overall agreement between theory and experiment (Fig. 2) is quite satisfactory. The agreement 
is less satisfactory for u = 32.5 meV and should be attributed to the breakdown of scaling since 
ujq > 7/2, as also evident in Fig. 3. Note that the ad hoc ansatz for f(x) commonly used 
(with A2 = 0) can be easily obtained assuming a Lorentzian dependence of Xq on A simple 



Figure 4. Left pannel: the temperature 
dependence of k 2 xq fo r several ( as a function of 
T. Right pannel: T-dependence of Xq relative to a 
Lorentzian C n . 



calculation yields arctan(Aiw/T) with A% ~ l/a, provided that k <C 1 but k 2 Xq ~ const (see 
Fig. 4). 



3. CP response: 

The advantage of the memory-function formalism is that the emergence of the CP at T < T g in 
the spin response can be as well treated within the same framework. One has to assume that 
unlike in a paramagnet the mode damping 7 q is not constant but may acquire an additional low 
frequncy contribution. In particular we can take M q ~ ij — 5 2 /(lo + iX), with T-dependent 8 
and A, which leads to XqO^) of the form used also to analyse experimental INS data for YBCO 
with x = 0.35 II]. 

For A — > but 5 2 /X S> 7 the modified M q leads to two distinct energy scales and hence to 
two contributions to spin dynamics, i.e., the CP part X q ( w ) an d the regular contribution x q ( w )> 
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valid for u < A and A < u -C 7, respectively. Thus, below T 5 new scales are set by T r = ^q/7 
and T c = (^/5 2 )A/x q , with x q ° = ??/^q and f2 2 = w 2 + (5 2 . If one assumes that 5 saturates 
at low T, as is manifest by saturation of r r |5j, r c oc A/x q becomes the smallest energy scale, 
resulting in a quasielastic peak of width T c . The saturation of r r , although at present unclear 
physically, is responsible for the transfer of spectral weight, since Cq ~ C q — irT/^Tr) [TO] , 
which is again consistent with experiment on x = 0.35 YBCO [5j. 



4. Conclusions 

The approach presented gives a consistent explanation of the ui/T scaling both in Xq( w ) as wen 
as in Xl( w )- It is based on two well established experimental facts: the overdamped nature of 
the response and the saturation of kt at low u and T. The appearance of the CP for T < T g is 
easily incorporated into the formalism via the (almost) singular oo- and T-dependent damping 
Mq(u)). However, the question as to the origin of CP remains to be settled. 
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